[Identification and expression analysis of whole gene family of Isatis indigotica 4-coumarate: CoA ligase].

The 4-coumarate: CoA ligase(4CL) is a key enzyme in the upstream pathway of phenylpropanoids such as flavonoids, soluble phenolic esters, lignans, and lignins in plants. In this study, 13 4CL family members of Arabidopsis thaliana were used as reference sequences to identify the 4CL gene family candidate members of Isatis indigotica from the reported I. indigotica genome. Further bioinformatics analysis and analysis of the expression pattern of 4CL genes and the accumulation pattern of flavonoids were carried out. Thirteen 4CL genes were obtained, named Ii4CL1-Ii4CL13, which were distributed on chromosomes 1, 2, 3, 4, and 6. The analysis of the gene structure and conserved structural domains revealed the intron number of I. indigotica 4CL genes was between 1 and 12 and the protein structural domains were highly conserved. Cis-acting element analysis showed that there were multiple response elements in the promoter sequence of I. indigotica 4CL gene family, and jasmonic acid had the largest number of reaction elements. The collinearity analysis showed that there was a close relationship between the 4CL gene family members of I. indigotica and A. thaliana. As revealed by qPCR results, the expression analysis of the 4CL gene family showed that 10 4CL genes had higher expression levels in the aboveground part of I. indigotica. The content assay of flavonoids in different parts of I. indigotica showed that flavonoids were mainly accumulated in the aboveground part of plants. This study provides a basis for further investigating the roles of the 4CL gene family involved in the biosynthesis of flavonoids in I. indigotica.

[Dead heart of pith-decayed Scutellariae Radix: a study based on multi-omics].

Dead heart is an important trait of pith-decayed Scutellariae Radix. The purpose of this study was to clarify the scientific connotation of the dead heart using multi-omics. Metabolomics and transcriptomics combined with multivariate statistical analysis such as principal component analysis(PCA) and partial least squares discriminant analysis(PLS-DA) were used to systematically compare the differences in chemical composition and gene expression among phloem, outer xylem and near-dead xylem of pith-decayed Scutella-riae Radix. The results revealed significant differences in the contents of flavonoid glycosides and aglycones among the three parts. Compared with phloem and outer xylem, near-dead xylem had markedly lowered content of flavonoid glycosides(including baicalin, norwogonin-7-O-ß-D-glucuronide, oroxylin A-7-O-ß-D-glucuronide, and wogonoside) while markedly increased content of aglycones(including 3,5,7,2',6'-pentahydroxy dihydroflavone, baicalin, wogonin, and oroxylin A). The differentially expressed genes were mainly concentrated in KEGG pathways such as phenylpropanoid metabolism, flavonoid biosynthesis, ABC transporter, and plant MAPK signal transduction pathway. This study systematically elucidated the material basis of the dead heart of pith-decayed Scutellariae Radix with multiple growing years. Specifically, the content of flavonoid aglycones was significantly increased in the near-dead xylem, and the gene expression of metabolic pathways such as flavonoid glycoside hydrolysis, interxylary cork development and programmed apoptosis was significantly up-regulated. This study provided a theoretical basis for guiding the high-quality production of pith-decayed Scutellariae Radix.

[Cloning and catalytic analysis of Isatis indigotica chalcone isomerase in vitro].

Chalcone isomerase is a key rate-limiting enzyme in the biosynthesis of flavonoids in higher plants, which determines the production of flavonoids in plants. In this study, RNA was extracted from different parts of Isatis indigotica and reverse-transcribed into cDNA. Specific primers with enzyme restriction sites were designed, and a chalcone isomerase gene was cloned from I. indigotica, named IiCHI. IiCHI was 756 bp in length, containing a complete open reading frame and encoding 251 amino acids. Homology analysis showed that IiCHI was closely related to CHI protein of Arabidopsis thaliana and had typical active sites of chalcone isomerase. Phylogenetic tree analysis showed that IiCHI was classified into type Ⅰ CHI clade. Recombinant prokaryotic expression vector pET28a-IiCHI was constructed and purified to obtain IiCHI recombinant protein. In vitro enzymatic analysis showed that the IiCHI protein could convert naringenin chalcone into naringenin, but could not catalyze the production of liquiritigenin by isoliquiritigenin. The results of real-time quantitative polymerase chain reaction(qPCR) showed that the expression level of IiCHI in the aboveground parts was higher than that in the underground parts and the expression level was the highest in the flowers of the aboveground parts, followed by leaves and stems, and no expression was observed in the roots and rhizomes of the underground parts. This study has confirmed the function of chalcone isomerase in I. indigotica and provided references for the biosynthesis of flavonoid components.

[Enzymatic characterization of lignan glucosyltransferase of Isatis indigotica].

The lignan glycosyltransferase UGT236(belonging to the UGT71 B family) from Isatis indigotica can catalyze the production of phloridzin from phloretin in vitro. UGT236 shares high identity with P2'GT from apple. In this study, the recombinant plasmid pET28 a-MBP-UGT236 was transferred into Escherichia coli Rosetta(DE3) cells and induced by isopropyl-ß-D-thiogalactoside(IPTG). The purified UGT236 protein was used for enzymatic characterization with phloretin as substrate. The results showed that UGT236 had the optimal reaction temperature of 40 ℃ and the optimal pH 8(Na_2HPO_4-NaH_2PO_4 system). The UGT236 activity was inhibited by Ni~(2+) and Al~(3+), enhanced by Fe~(2+), Co~(2+), and Mn~(2+), and did not affected by Mg~(2+), Ca~(2+), Li~+, Na~+, or K~+. The K_m, K_(cat), and K_(cat)/K_m of phloretin were 61.03 µmol·L~(-1), 0.01 s~(-1), and 157.11 mol~(-1)·s~(-1)·L, and those of UDPG were 183.6 µmol·L~(-1), 0.01 s~(-1), and 51.91 mol~(-1)·s~(-1)·L, respectively. The possible active sites were predicted by homologous modeling and molecular docking. By mutagenisis and catalytic activity detection, three key active sites, Glu391, His15, and Thr141, were identified, while Phe146 was related to product diversity. In summary, we found that the lignan glycosyltransferase UGT236 from I.indigotica could catalyze the reaction of phloretin into phloridzin. Several key amino acid residues were identified by structure prediction, molecular docking, and site-mutagenesis, which provided a basis for studying the specificity and diversity of phloretin glycoside products. This study can provide a reference for artificially producing glycosyltransferase elements with high efficiency and specific catalysis.

[Cloning and functional characterization of lignan glycosyltransferase gene IiUGT349 in Isatis indigotica].

Based on the transcriptome data of Isatis indigotica, a total of 110 putative glycosytransferases were identified. Through prokaryotic expression and enzymic activity assay in vitro, a novel lignan glycosyltransferase gene was screened out and named IiUGT349, which catalyzed lariciresinol into lariciresinol-4-O-ß-D-glucoside and lariciresinol-4'-O-ß-D-glucoside. Bioinformatics analysis suggested that IiUGT349 contained an open reading frame(ORF) of 1 401 bp encoding a protein of 467 amino acids. A protein analysis indicated that IiUGT349 have a predecited molecular weight of 52.77 kDa and pI of 5.96. Phylogenetic analysis showed that IiUGT349 belonging to UGT90 family shared low amino acid sequence identity with the reported lignan glycosyltransferases, which may represent a novel type of lignan glycosyltransferases. Quantitative real-time PCR(qRT-PCR) analysis showed that IiUGT349 was expressed in roots, stems, young leaves and leaves, with the highest expression level in stems. Further biochemical analysis showed that the optimal reaction time of IiUGT349 recombinant protein was 12 h and the optimal temperature was 45 ℃. Subcellular localization demonstrated that IiUGT349 was located in the cytoplasm and nucleus of plants. In this study, a new glucosyltransferase gene IiUGT349 from I. indigotica belonging to the UGT90 family was cloned, which laid a foundation to further investigate its' function and elucidate the lignan glycosides biosynthesis pathway and plays an important role for great significance for the synthetic biology of active lignan glycosides.

Real-time toxicity prediction of Aconitum stewing system using extractive electrospray ionization mass spectrometry.

Due to numerous obstacles such as complex matrices, real-time monitoring of complex reaction systems (e.g., medicinal herb stewing system) has always been a challenge though great values for safe and rational use of drugs. Herein, facilitated by the potential ability on the tolerance of complex matrices of extractive electrospray ionization mass spectrometry, a device was established to realize continuous sampling and real-time quantitative analysis of herb stewing system for the first time. A complete analytical strategy, including data acquisition, data mining, and data evaluation was proposed and implemented with overcoming the usual difficulties in real-time mass spectrometry quantification. The complex Fuzi (the lateral root of Aconitum)-meat stewing systems were real-timely monitored in 150 min by qualitative and quantitative analysis of the nine key alkaloids accurately. The results showed that the strategy worked perfectly and the toxicity of the systems were evaluated and predicated accordingly. Stewing with trotters effectively accelerated the detoxification of Fuzi soup and reduced the overall toxicity to 68%, which was recommended to be used practically for treating rheumatic arthritis and enhancing immunity. The established strategy was versatile, simple, and accurate, which would have a wide application prospect in real-time analysis and evaluation of various complex reaction systems.

Glucosyltransferase Capable of Catalyzing the Last Step in Neoandrographolide Biosynthesis.

ApUGT, a diterpene glycosyltransferase from Andrographis paniculata, could transfer a glucose to the C-19 hydroxyl moiety of andrograpanin to form neoandrographolide. This glycosyltransferase has a broad substrate scope, and it can glycosylate 26 natural and unnatural compounds of different structural types. This study provides a basis for exploring the glycosylation mechanism of ent-labdane-type diterpenes and plays an important role in diversifying the structures used in drug discovery.

Chemical profiling in Moutan Cortex after sulfuring and desulfuring processes reveals further insights into the quality control of TCMs by nontargeted metabolomic analysis.

As a traditional processing method, sulfuring has been used in the processing of many traditional Chinese medicines (TCMs). Desulfuring, which has emerged in recent years, is a new method applied to sulfured herbs so they can comply with regulations regarding residual SO2. Due to the chemical transformations and the residual SO2 in the herbs, both sulfuring and desulfuring have negative effects on the safety and therapeutic effects of TCMs, and Moutan Cortex is one of the TCMs most susceptible to these effects. Here, a new strategy was developed to differentiate normal, sulfured and desulfured Moutan Cortex, and the transformations of compounds in sulfuring and desulfuring processes were analyzed using ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC/Q-TOF-MSE) method based on metabolomic analysis. Our findings were as follows: (1) a total of 119 compounds were identified or tentatively identified, including 9 compounds that are being reported for the first time as natural products; (2) 15 sulfocompounds were generated during the sulfuring process; (3) these sulfocompounds could not be converted back into their corresponding glycosides by the desulfuring process, and the desulfuring decreased the residual SO2,while also removing some soluble compounds in the sulfured Moutan Cortex; and (4) 28 compounds were screened and tentatively identified as markers for distinguishing normal, sulfured and desulfured Moutan Cortex. Our findings provide a new practical strategy for evaluating how sulfuring and desulfuring affect the quality of TCMs.

Structural characterization and discrimination of the Paris polyphylla var. yunnanensis and Paris vietnamensis based on metabolite profiling analysis.

This study aimed to distinguish the rhizomes of Paris polyphylla var. yunnanensis (Franch) Hand Mazz (PPY) and Paris veitnamensis (Takht.) H. Li (PV) using metabolomics-based ultra high-performance liquid chromatography coupled with quadrupole time-of-fligh mass spectrometry (UHPLC/Q-TOF MS). First, the UHPLC/Q-TOF MS approach was optimized for metabolite profiling. Then, the MS data were processed using UNIFI™ combined with an in-house library to automatically characterize the metabolites. Based on the exact mass information, the fragmentation characteristics, and the retention time of compounds, and the fragmentation mechanism and retention behavior of steroidal glycosides in the references, the structures identified by UNIFI were further verified. Overall, 146 metabolites, including 42 potential new compounds, were identified or tentatively identified. Pattern recognition analysis of the PPY and PV MS data revealed that they were clearly separated, and 15 potential biomarkers for differentiating between them were selected. These biomarkers were subsequently used to successfully predict the genus of PPY and PV samples. These results indicated that metabolite profiling by UHPLC/Q-TOF MS is an effective, robust approach for determining the characteristic biomarkers that differentiate between TCM species with multiple botanical origins.

[Mitogen-activated protein kinase genes of Artemisia annua and their expression analysis under Cadmium stress].

In order to study Artemisia annua under cadmium stress, whether there are corresponding MAPK genes involved in transduction of the cadmium signal. 17 AaMAPK genes, named AaMAPK1-AaMAPK17 repectively, were finally obtained by using Trinity method for de novo assembly of transcripts from SRA database and BLAST search against AtMAPK genes and determing conserved domain using a series of bioinformatics tools. There exist 16 MAPK genes contained T[D/E]Y conserved domains among the obtained genes. The expressions of these genes were analyzed by Real-time PCR under cadmium stress. The results showed that the expressions level of AaMAPK3 and AaMAPK10 were down-regulated and MAPK7, MAPK9 and MAPK12 were up-regulated. These indicated that there exist corresponding MAPK genes involved in transduction of the cadmium signal.

[Status and future of natural resource for Chinese materia medica].

For thousands of years, the natural resource for Chinese materiamedica has been the foundation of the traditional Chinese medicine industry, which provides abundant medicine for human. In recent years, increasing demands and irrational exploitation led to a lot of problems such as rapid decrease of traditional Chinese herbs reserves, low quality of medicine and dismishing traditional cultures. These restricted the development of the traditional Chinese medicine. To solve these problems, scientists have done much work on investigating traditional Chinese medicine resources, exploring the metabolic pathway of bioactive ingredients, cultivating new varieties, and carrying out synthetic biology. These studies provided a theoretical basis for sustainable utilizationand future developmentof traditional Chinese medicine resources.

The Biosynthetic Pathways of Tanshinones and Phenolic Acids in Salvia miltiorrhiza.

Secondary metabolites from plants play key roles in human medicine and chemical industries. Due to limited accumulation of secondary metabolites in plants and their important roles, characterization of key enzymes involved in biosynthetic pathway will enable metabolic engineering or synthetic biology to improve or produce the compounds in plants or microorganisms, which provides an alternative for production of these valuable compounds. Salvia miltiorrhiza, containing tanshinones and phenolic acids as its active compounds, has been widely used for the treatment of cardiovascular and cerebrovascular diseases. The biosynthetic analysis of secondary metabolites in S. miltiorrhiza has made great progress due to the successful genetic transformation system, simplified hairy roots system, and high-throughput sequencing. The cloned genes in S. miltiorrhiza had provided references for functional characterization of the post-modification steps involved in biosynthesis of tanshinones and phenolic acids, and further utilization of these steps in metabolic engineering. The strategies used in these studies could provide solid foundation for elucidation of biosynthetic pathways of diterpenoids and phenolic acids in other species. The present review systematically summarizes recent advances in biosynthetic pathway analysis of tanshinones and phenolic acids as well as synthetic biology and metabolic engineering applications of the rate-limiting genes involved in the secondary metabolism in S. miltiorrhiza.

[Cloning and expression analysis of a hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferasegene(IiHCT) from Isatis indigotica].

Based on the transcriptome data, we cloned the open reading frame of IiHCT gene from Isatis indigotica, and then performed bioinformatic analysis of the sequence. Further, we detected expression pattern in specific organs and hairy roots treated methyl jasmonate( MeJA) by RT-PCR. The IiHCT gene contains a 1 290 bp open reading frame( ORF) encoding a polypeptide of 430 amino acids. The predicted isoelectric point( pI) was 5.7, a calculated molecular weight was about 47.68 kDa. IiHCT was mainly expressed in stem and undetectable in young root, leaf and flower bud. After the treatment of MeJA, the relative expression level of IiHCT increased rapidly. The expression level of IiHCT was the highest at 4 h and maintained two fold to control during 24 h. In this study, cloning of IiHCT laid the foundation for illustrating the biosynthesis mechanism of phenylpropanoids in I. indigotica.

[Stability analysis of reference gene based on real-time PCR in Artemisia annua under cadmium treatment].

In this study, Actin, 18S rRNA, PAL, GAPDH and CPR of Artemisia annua were selected as candidate reference genes, and their gene-specific primers for real-time PCR were designed, then geNorm, NormFinder, BestKeeper, Delta CT and RefFinder were used to evaluate their expression stability in the leaves of A. annua under treatment of different concentrations of Cd, with the purpose of finding a reliable reference gene to ensure the reliability of gene-expression analysis. The results showed that there were some significant differences among the candidate reference genes under different treatments and the order of expression stability of candidate reference gene was Actin > 18S rRNA > PAL > GAPDH > CPR. These results suggested that Actin, 18S rRNA and PAL could be used as ideal reference genes of gene expression analysis in A. annua and multiple internal control genes were adopted for results calibration. In addition, differences in expression stability of candidate reference genes in the leaves of A. annua under the same concentrations of Cd were observed, which suggested that the screening of candidate reference genes was needed even under the same treatment. To our best knowledge, this study for the first time provided the ideal reference genes under Cd treatment in the leaves of A. annua and offered reference for the gene expression analysis of A. annua under other conditions.

[In silico cloning and bioinformatics analysis of an AP2/EFR family gene from Arnebia euchroma].

A cDNA sequence of Arnebia euchroma AP2/ERF named AeAP2/ERF was cloned by in silico cloning in this study, using ACX71873 sequence from Lithospermum erythrorhizon as the probe sequence. Some characters of the AP2/ERF gene and encoded protein sequences were predicted and analyzed by the bioinformatics methods, including general physical and chemical properties, hydrophobieity, signal peptide, secondary structure, localization sites in cells. Results showed that the 876 bp long gene included a 1 077 bp ORF and encoding 205 amino acid. The AeAP2/ERF protein had no signal peptide, it was a hydrophilic proteins located in nucleus. The function of the AP2/ERF protein was mainly involved with metabolism controlling and signal transduction.